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Abstract 

The origin and evolution of Lithium-Beryllium-Boron is a crossing point be- 
tween different astrophysical fields : optical and gamma spectroscopy, non ther- 
mal nucleosynthesis, Big Bang and stellar nucleosynthesis and finally galactic 
evolution. We describe the production and the evolution of Lithium-Beryllium- 
Boron from Big Bang up to now through the interaction of the Standard Galactic 
Cosmic Rays with the interstellar medium, supernova neutrino spallation and a 
low energy component related to supernova explosions in galactic superbubbles. 
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1 Introduction 




Figure 1: Big Bang Nucleosynthesis of Lithium, BeryUium and Boron vs the 
photon over baryon ratio. N/H represents the abundance by number of the 



different nuclear species 92 



Light element nucleosynthesis is an important chapter of nuclear astro- 
physics. Specifically, the rare and fragile light nuclei, Lithium, Beryllium and 
Boron (LiBeB) are not generated in the normal course of stellar nucleosynthesis 
(except ^Li, in the galactic disk) and are, in fact, destroyed in stellar interiors. 
This characteristic is reflected in the low abundance of these simple species. 

A glance to the abundance curve suffices to capture the essence of the 
problem: a gap separates He and C. At the bottom of this precipice rests the 
trio Li-Be-B. They are characterized by the simplicity of their nuclear structure 
(6 to 11 nucleons) and their scarcity in the Solar System and in stars. 

Indeed, they are rare because they are fragile and, apparently a selection 
principle at the nuclear level has operated in nature. Due to the fact that 
nuclei with mass 5 and 8 are unstable, the Big-Bang nucleosynthesis (BBN) 
has stopped at A = 7, and primordial thermonuclear fusion has been unable to 
proceed efhciently beyond lithium. Figure 1 represents the calculated abundance 
by number of the light elements as a function of the baryon over photon ratio 
produced in the Big Bang. 

The Big Bang production of ®Li is dominated by the D (a, 7) ®Li reaction 
[ p8| , , , . No direct measurement of the cross section of this reaction 
has been performed below 1 MeV. However, the Coulomb breakup technique 
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[ p^ provides an indirect estimate which is in qualitative agreement with the 
theoretical extrapolation at low energy of Mohr et al. (1994) Recently, 
the European Collaboration between nuclear physicits and astrophysicits led by 
Marcel Arnould (NACRE : European Astrophysical Compilation of Reaction 
Rates) has delivered a consistent compilation of thermonuclear reaction rates 
of astrophysical interest, among them is the D (0,7) ^Li reaction [Q. They 
conclude that the reaction rate based on the Mohr et al. S factor is the 
most relevant. This rate is similar to that of Caughlan and Fowler (1988) 
in the temperature range of cosmological interest. The two estimates agreee 
to within a few percents. NoUet et al. (1997) |]l^ have considered all the pub- 
lished evaluations of this reaction rate. However, most of the extrapolations at 
low temperature depart considerably from the Kiener et al. 1991 [^2] estimate 
except that of ||l^. Following the recommendation of Kiener (1998, private com- 
munication) and we have adopted the Mohr estimate. Note that the upper 
limit given by ^ is much higher. This upper limit is indeed related to the bad 
sensitivity of the detector used (Kiener, private communication). The and 
^^B abundances are calculated with updated reactions rates, including the new 
^°B(p, a)''Be reaction ( adopted from (l5)). Clearly, the calculated primordial 
Be and B abundances are negligibly small, compared to ^Li. The standard BBN 
is hoplessly ineffective in generating ^Li, ^Be, ^°B, ^^B. iQ, |11 ). 

Thus, stars are necessary to pursue the nuclear evolution bridging the gap 
between ''He and ^^C much later, through nuclear fusion. 

Up to recently, the most plausible formation agents of LiBeB were thought 
to be Galactic Cosmic Rays (GCRs) interacting with interstellar CNO. Other 
possible origins have been also identified: primordial and stellar (''Li) and su- 
pernova neutrino spallation (for ''Li and ^'B). In contrast, ^Li, ^Be and "^B are 
pure spallative products. Be is very precious to astrophysics since it is monoiso- 
topic, thus the elemental measurement is also an isotopic measurement. On the 
other hand, Li and B have two stable isotopes and an isotopic measurement is 
necessary to separate 6 and 7, 10 and 11, which is very difficult, specifically in 
stars. 

^Li presents a special interest since fortunatly, the ^Li/^Li ratio has been 
measured recently in a few halo stars, offering a new constraint on the early 
galactic evolution of light elements. Optical measurements of the beryllium and 
boron abundances in halo stars have been achieved by the 10 meter KECK 
telescope and the Hubble Space Telescope. These observations indicate a quasi 
linear correlation between Be and B vs Fe, at least at low metallicity (mass 
fraction of elements heavier than helium), contradictory at first sight, to a dom- 
inating GCR origin of the light elements which predicts a quadratic relationship 
(see appendix, section 10). As a consequence, the theory of the origin and evo- 
lution of LiBeB nuclei has to be reassessed. Aside GCRs, which are thought 
to be accelerated in the general interstellar medium (ISM) and which create 
LiBcB through the break up of interstellar CNO by their fast protons and al- 
phas, Wolf-Rayet stars (WR) and core collapse supernovae (SNII) grouped in 
superbubbles could produce copious amounts of light elements via the fragmen- 
tation in flight of rapid carbon and oxygen nuclei (called hereafter low energy 
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component, LEG) colliding with H and He in the ISM. In this case, LiBeB would 
be produced independently of the interstellar medium chemical composition and 
thus a primary origin is expected (see appendix, section 10). These different 
formation processes are discussed in the framework of a galactic evolutionary 
model. More spectroscopic observations (specifically of O, Fe, Li, Be, B) in halo 
stars are required for a better understanding of the relative contribution of the 
various mechanisms. Future tests on the injection and acceleration of nuclei by 
supernovae and Wolf Rayet relying on high energy astronomy will be invoked 
in the perspective of X ray astronomy and the European INTEGRAL satellite. 

Thus, light element research impacts several important astrophysical prob- 
lems, specifically, the origin and the evolution of cosmic rays, galactic chemical 
evolution, X and gamma-ray astronomies and indirectly BBN through Li. 

The questions posed are: What are the main agents of nonthermal nucle- 
osynthesis and what are their relative contributions at different galactic epochs? 
What is the origin of the present epoch GGRs, are they accelerated out of the 
average ISM or out of supernova ejecta, and what was their past composition? 
Was there a population of low energy nuclei, specifically G and O, responsible 
for primary LiBeB production? 

This review is dedicated to David Schramm, our exemplary colleague and 
friend, who has largely contributed to the development of this astrophysical 
field. 

In section 2, we decline the story of the subject matter, in section 3 the 
production ratios of the light isotopes are presented. Sections 4 and 5 are 
devoted to standard Galactic Gosmic Rays, the traditional agents of creation of 
LiBeB. Section 6 introduces and describes a new spallative component made of 
low energy nuclei, possibly distinct from GGRs. Section 7 briefly summarises 
neutrino spallation, section 8 integrates all these processes in the framework 
of a galactic evolutionary model. Finally, section 9 concludes and presents the 
prospects of LiBeB research. 

2 LiBeB story 

A general trend in nature is that complex nuclei are not proliferating: the 
abundance of the elements versus the mass number draws a globally decreasing 
curve. In the whole nuclear realm, LiBeB are exceptional since they are both 
simple and rare. Typically, in the Solar System, Li/H = 2. 10~^, B/H = 7. 
10"^°, Be/H = 2.5 10"" 0, whereas G/H = 3.5 lO"* and 0/H = 8.5 lO"*. 

The local isotopic ratios, measured in meteorites are known with excellent 
precision. "B/^OB = 4 [||, ^Li/^Li = 12.5 §. 

Stellar nucleosynthesis, quiescent or explosive, forge the whole variety of 
nuclei from G to U but LiBeB nuclei are destroyed in the interior of stars, 
except ^Li which is produced in AGB and novae. The destruction temperatures 
are 2, 2.5, 3.5, 5.3 and 5 millions of degrees for ^Li, ^Li, ^Be, ^°B and ^^B 
respectively. It is worth noting that ^Li and ^^B could be produced by neutrino 
spallation in helium and carbon shells of core collapse supernovae, respectively 
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[ p6| , ; however, this mechanism is particularly uncertain depending strongly 
on the neutrino energy distribution. 

It is clear that another source is necessary to generate at least ^Li, ^Be, ^"B 
and this is a non thermal mechanism, namely the break up of heavier species 
(CNO, mainly) by energetic collisions, also called spallation. 

The LiBeB story has been rich and moving. The genesis of LiBeB was so 
obscure to Burbidge et al (1957) ||] that they called X the process leading to 
their production. Then came Hubert Reeves and his colleagues. 

A very active and talented group of nuclear physicists, led by Rene Bernas: 
Eli Gradstzjan, Robert Klapish, Marcelle Epherre, Francoise Yiou and Grant 
Raisbeck used the mass spectroscopy techniques in order to determine the spal- 
lation cross sections induced by 155 MeV protons on C and O. H. Reeves and 
E. Schatzman were lucky enough to collaborate with this group. The paper 
Bernas et al. (1967) jsj based on a better knowledge of the spallation nuclear 
reaction cross sections has been determining to the understanding of the genesis 
of LiBeB. 

The second happy circumstance was the encounter between Hubert Reeves 
and Bernard Peters (one of the discoverers of the heavy Cosmic Rays). B. 
Peters draw the attention of H. Reeves on the fact that the LiBeB/CNO ratio is 
about 10~^ in the standard abundances (Solar System) and 0.25 in the GCRs. 
H. Reeves proposed at once the scenario concerning the LiBeB nucleosynthesis 
[ [78| : contrary to most of the nuclear species LiBeB are formed by the " spallative 
" encounter between the energetic GCR particles and the interstellar medium 
(ISM). 

In a seminal work, Meneguzzi, Audouze and Reeves (1971) (MAR) gave 
the complete calculations. Considering the fast p, a in the GCRs interacting 
with CNO in the ISM, they were able to make quantitative estimates of the 
LiBeB production on the basis of cross section measurements notably made in 
Orsay 0. 

However, this study, relying on the local and present observations (LiBeB 
and CNO abundances, cosmic ray flux and spectrum) was based on an ex- 
trapolation over the whole galactic lifetime assuming that all the parameters 
are constant. This result accounted fairly well for the cumulated light element 
abundances but obviously not for their evolution which, at that time, was un- 
known. The pertinence of the MAR proposal is illuminated by the simple and 
beautiful fact that the hierarchy of the abundances ^^B > ^°B > ^Be is reflected 
in the cross sections |Q (see figure 4) . This is another proof that nature follows 
the rules of nuclear physics. ®Li, ^Be and ^°B were nicely explained but prob- 
lems were encountered with ^Li and ^^B. The calculated ^Li/^Li ratio was 1.2 
against 12.5 in meteorites. Stellar sources of ^Li appeared necessary. The esti- 
mated ^^B/^"B ratio was 2.5 instead of 4 in meteorites. An ad-hoc hypothesis 
drawing on unobservable low energy protons operating through the ^^N(p,x)^^B 
reaction was advocated ||7^ . 

New measurements of Be/H and B/H from KECK and HST, together with 



[Fe/H] [118|, |110| g |105|, |119 , [107|, |109| in very low metallicity halo stars 



came to set strong constraints on the origin and evolution of light isotopes. 
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Figure 2: Lithium (essentially ""Li) and Beryllium abundances as a function of 



the iron abundance of stars, kindly supplied by Martin Lemoine. |99, m, 111 



105, 110, 119| 



The evolution of BeB was suddenly uncovered over about 10 Gyr, taking 
[Fe/H] as an evolutionary index. A compilation of Be and B data is presented 
in figures 2 and 6. The most striking point is that log(Bc/H) and log(B/H) are 
both quasi proportional to [Fe/H] (this notation means the log of the number 
abundance normalized to its solar value), at least up to [Fe/H] = -1 and that the 
B/Be ratio lies in the range 10 - 30 ]107| . Note however, two discrepant points 
in the boron diagram (figure 6) at the lowest [Fe/H]. This is mainly due to the 
huge NLTE correction on B data ]112| ] , [Q that increases the departure from a 
straight line. It is important to take a carefuU look to this delicate correction. 
Moreover, NLTE corrections on Fe have also to be considered |88 . 

This linearity came as a surprise since a quadratic relation was expected 
from the GCR mechanism (see section 4). It was a strong indication that the 
standard GCRs are not the main producers of LiBeB in the early Galaxy. A 
new mechanism of primary nature was required to reproduce these observations: 
it has been proposed that low energy CO nuclei produced and accelerated by 
massive stars (WR and SN II) fragment on H and He at rest in the ISM. This 
low energy component (LEG) has the advantage of coproducing Be and B in 
good agreement with the ratio observed in Pop 11 stars (figure 6, |3^, [Q, | ]93{ , 

[0, II.)- . . . . 

A primary origin , in this language, means a production rate independent of 

the interstellar metallicity (Z: mass fraction of elements heavier than helium). In 
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this case, the cumulated abundance of a given light isotope L is approximately 
proportional to Z. At variance, standard GCRs offer a secondary mechanism 
because it should depend both on the CNO abundance of the ISM at a given 
time and on the intensity of cosmic ray flux, itself assumed to be proportional 
to the SN II rate (see appendix, section 10). 

Indeed, the Be-Fe and B-Fe correlations taken at face value reveal a contra- 
diction between standard GCR theory and observation. But, since oxygen is 
the main progenitor of BeB, the apparent linear relation between BeB and Fe 



could be misleading if O was not strictly proportional to Fe , . Thus the 



pure primary origin of BeB in the early Galaxy could be questionned , |^ , 
[ p5| . However, the oxygen measurements themselves are confronted to many 
difficulties 114 |, [2^ , fs^ . On the theoretical side, the sit uati on is not better. 
The [a/Fe] vs [Fe/H] where a = Mg, Si, Ca, S, Ti |l0|l, |l2| show a plateau 
from [Fe/H] = -4 to -1. On nucleosynthetic grounds, it would be surprising that 
oxygen does not follow the Si a nd Ca trends. Moreover, using the published 
nucleosynthetic yields |9^, |l2l| it is impossible through galactic evolutionary 
models to fit the log(0/H) vs [Fe/H] relation of Israehan et al. (1998) [^ and 
Boesgaard et al. (1999) |Q since the required oxygen yields are unrealistic. Thus 
the subject is controversial. 

Concerning lithium, a compilation of the data is shown in Figure 2 |]113| and 



[117|. The Spite plateau extends up to [Fe/H] = -1.3. Beyond, Li/H is strongly 
increasing until its solar value of 2. 10~^. 

A stringent constraint on any theory of Li evolution is avoiding to cross the 
Spite's plateau below [Fe/H] = -1. Accordingly, the Li/Be production ratio 
should be less than about 100. 

Recent measurements of ^Li have been made successfully in two halo stars, 
HD84937 and BD +26 3578 at about [Fe/H] = -2.3 Q |l0| pO| ], pT[ , 

yielding ^Li/^Li about 0.05. The great interest of ^Li, (as shown by |102| ) 
besides of being an indicator of stellar destruction [lOl], ||T^, [Q, ]123| ] is to 
represent a pure spallation product as ^Be and participate to constrain the 
global LiBeB evolution. Moreover, the ^Li/^Be ratio (clearly non solar in two 
halo stars since its amounts to 20 - 80, |28|] ) bears information on ancient non 
thermal nuclei. 

To summarize, we can give six observational constraints on LiBeB evolution: 

1. Be and B proportional to Fe 

2. Li/Bc < 100 up to [Fc/H] = -1 

3. B/Bc = 10-30 

4. "B/i"B 4 at solar birth 

5. ■^Li/^Li = 12.5 at solar birth 

6. SLi/'^Li = 0.05 and ^Li/^Bc = 20 to 80 (to be confirmed) at [Fe/H] about 
-2.3 

We recall that the observational O - Fe relation is central to the interpreta- 
tion since specifically the production of Be is directly related to O and not to 
Fe. 
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3 Production ratios of LiBeB isotopes 



3.1 Cross sections 

The physics of interest for us is that of the heavy ions accelerators. It implies 
a beam of fast nuclei (injected and subsequently accelerated), a target and an 
interaction between both. 

The result of a nuclear collision depends on i) the composition of the beam, ii) 
the relative velocity of the projectile and the target nuclei, and thus the energy 
spectrum of the beam and iii) the target composition. All these dependencies 
are quantified by the cross section. 

For a moncncrgctic beam, the number of favourable events per second or 
reaction rate is given by = riirijav, or, equivalently, by the product riia.cj) 
where rii is the number density of the target and (f) the flux of the projectiles. For 
an energy distribution, the reaction rate is averaged on the velocity distribution 
function </'(v). In most of the astrophysical situations, the energy distribution 
of the projectile is violently non thermal. The spectra that can be parametrized 
as power laws of the form 

N(E)dE = kE-'-'dE 

are customary. Stellar nucleosynthesis implies low energies (1 keV-100 keV) 

and high densities (10^ — 10^ g.cm~^). By contrast, spallativc nucleosynthe- 
sis implies high energies (1 MeV-100 GeV) and very low densities (1- 10^ 
atoms. cm~^). Expressed in (MeV/n) unit, the cross sections of direct reactions 
(1 + 2 — > 3) and reverse (2 + 1 — > 3) where beam and target have been 
exchanged, are identical due to the velocity conservation in energetic collisions. 

In the thermonuclear context, the energy of the nuclei is below the Coulomb 
barrier, and the cross section is dominated by the Coulomb penetration factor. 
Cross sections are factorised as follows: 

cr{E) = {l/E).S{E).exp{-2'Kri) 

where 77 = ZiZ2/hv 

The first term (square of the associated wavelength) is of geometric nature, 
the second, called astrophysical factor, is related to the internal structure of the 
target nucleus and the third, the Coulomb penetration factor, is by far the most 
influential. 

In general, the production cross sections, and thus the production ratios of 

the various isotopes are very different, which explains the abundance disparities 
at the end of the Big Bang nucleosynthesis. In the high energy non thermal 
context, on the contrary, the penetration factor tends to 1. This, added to the 
nuclear similarity of light nuclei leads to spallative production rates of ^Li,^Li, 
^Be and ^°B, ^^B relatively similar. 

The excitation functions (variations of the cross sections as a function of 
energy/n) are presented in figure 3. 

The a + a reaction, leading to the synthesis of ^Li and ^Li plays a peculiar 
role in astrophysics since it is the only reaction that implies nuclei of Big Bang 
origin. It is fertile even at zero mctallicity and thus it is a source of ®Li and ''Li 
in the early Galaxy. Its excitation function is particular due to its low threshold 
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E(MeV/nucl) 

Figure 3: LiBeB production cross sections as a function of the energy per nuclcon 
updated by From top to bottom and left to right, reactions on C, N, and 
O induced by protons, a + a ; same reactions but with alpha particles. These 
graphs have been kindly supplied by Reuven Ramaty. 

and its decline at high energy. For more details on the nuclear physics see Reeves 
(1994), 0) 

3.2 Input parameters 

Thus, four parameters are influential to the spallative production of light el- 
ements: the reaction cross sections, the energy spectrum of fast nuclei, the 
composition of the beam and the composition of the target. 

Cross sections are well measured and have been updated recently by 

IE- 

The adopted spectra are generally of two kinds: 

1. GCR: N(E)dE = kE~'^ ''dE above a few GeV/n with a flatening below 
(e.g. H and @). 

2. LEG : Shock wave acceleration with a cut at Eo of the form N(E) dE = 
kE-^-^ exp(-E/Eo) dE [|||, propagated in the ISM. 

The source composition of GCR is well determined ]108| . It is p and a rich 
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(H/0 ~ 200, He/0 = 20) contrary to the possible source composition of the 
LEG (see section 6). The most obvious contributors to LEG are supernovae, 
Wolf-Rayet and mass loosing stars Q , Q , Q . It is worth noting that in the 
early Galaxy, supernovae play a leading role since at very low metallicity stellar 
winds are unsignificant. 

Table 1 shows a sample of compositions used by different authors: solar 
syste m (SS) for comparison from co smic ray source (GRS) (|Q from 
[I15|), wind of massive stars (W40) (I^Sffrom |ll6| ), composition of grain prod- 
ucts (GR) (|7| and 40 Mo supernova at Z = 10"^ Zo from , 35 Mo 
supernova of solar metallicity (j?^ from |124|). The two supernovae, though at 
different metallicities, (SN40 at low metallicity and SN35 at solar metallicity), 
give similar yields due to the fact that metallicity dependent mass loss has not 
been taken into account in the stellar models. Resulting elemental and isotopic 
ratios (B/Be^B/^^B, ^Li/^Be) for different compositions and Eo can be found 
in |7| and (l2|. 

Note that N is unsignificant and that Type II supernovae are 0-rich whereas 
the winds of massive stars are G and He-rich. These abundance differences 
are important since the highest ^^B/^^B ratios are produced by G-rich beams 
through^^G(p,x)-'^^B and the highest ^Li/^Be ratios are produced by He and O 
rich compositions. 



o 



101 

8 
7 
6 
5 



10" 
100 





(a) 




ChR 






GR 








Strong Shock 


GCR 


Ambient: SS 





101 



102 




EoCMeV/nucleon] 



10* 

S-i 

QJ 

CO 102 

ID 

o 

^101 
100 



10-1 









lOOM.ofB : 
Strong Shock ^ 
Ambient: SS : 










1 axio'SN' 





100 



m 

CD 



sl02 r 



'101 , 



100 

100 



IQi 



102 









strong Shock 








Ambient: SS 








--%5"--.. (d) ; 














ffC 


AG 






SN60 


— =-^-=i.,:,,^^^^^_GCRi 



101 



102 



Eo[MeV/nucleon] 



Figure 4: Production ratios of "B/1"B (a), B/Be (c) and ^Li/^Be (d) and 
energetics (b) as a function of cut-off energy Eo in the strong shock case, kindly 
supplied by Reuven Ramaty |Q . 

Eigure 4 shows a sample of production ratios together with the associated 
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energetics for different source compositions as a function of the cut-off energy 
in the case of a strong shock spectrum. The values of the ratios at high Eo 
converge toward the GCR points. The data (ChR, FOS, AG) are from || Q 
and [Q. The source compositions are following. CRS: cosmic ray source, SS: 
solar system, WC: Wolf-Rayet, SN35, SN60: supernovae with 35 and 60 Mo 
progenitors, GR: grains. Figure (4.b) represents the energy injected under the 
form of fast nuclei necessary to produce 100 Mo of boron (which is in order of 
magnitude the mass of boron in the whole Galaxy). The target composition is 
that of the present ISM taken as the solar one. The horizontal bar indicates 
the total energy released by 2.10'' SN, which corresponds roughly to the total 
number of SN having exploded since the birth of the Galaxy. 



Table 1 : Source Composition 



Element 


SS 


CRS 


W40 


GR 


SN40(low Z) 


SN35(Zo) 


H 


1200 


220 


80 


2 


37 


27 


He 


120 


22 


25 





8.8 


7.6 


C 


0.47 


0.87 


1.6 


0.3 


0.09 


0.08 


N 


0.13 


0.04 




0.03 









1 


1 


1 


1 


1 


1 



The fourth parameter, i.e., the composition of the target (ISM) varies from 



the birth of the Galaxy up to now. The extensive study of |7^ and |122] shows 
that there are only slight differences in the results when the ISM metallicity is 
varied between Z=0 (early galaxy) and Z=0.02 (at present). 

4 Standard Galactic Cosmic Rays (GCRs): ac- 
celeration and propagation 

4.1 High energy cosmic rays 

Galactic cosmic rays represent the only sample of matter originating from be- 
yond the Solar System. They are constituted by bare energetic nuclei stripped 
from their electrons. Their energy density (about leVcm~^ |Q| similar to that 
of stellar light and that of the galactic magnetic field), indicate that they are an 
important component in the dynamics of the Galaxy. A key point for our pur- 
pose is that, as said above, GCRs are exceptionnally LiBeB rich (LiBeB/CNO 
= 0.25) compared to the local (Solar System) matter {LiBeB/CNO = 10^*^). 

The energy spectrum of the GCRs can be described by a power law above 
a few GeV/n (section 3.2). Below IGeV/n, the nuclei are repelled by the solar 
wind and it is difficult to deduce from the observations near earth, the inter- 
stellar spectrum. 

Diffuse shock acceleration is now the leading paradigm in cosmic ray physics 
(for a review see [|o|, @). In this mechanism, each passage of the shock 
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front produces an energy increment. Diffusion allows confinement in the shock 
region, and after many crossings of the shock (both ways), the particles acquire 
sufficient energy (gyroradius) to leak out. 

The energy required to sustain the GCR energy density can be supplied by 
supernovae Q. Thus SN are not necessary the sources of GCR nuclei, but 
their accceleration agents. Recurrently models assuming that they are both the 
reservoir and the accelerating engine have been proposed. The last version of 
this idea |Q, though attractive has been criticized. This does not mean that 
fast nuclei originating from supernovae do not exist at all. If they exist they are 
different from GCR and should be called differently (see section 6). 

In the simplest (linear) shock acceleration theory, the resulting momentum 
distribution function f(p) is a power law of the form p~3r/(r-i)^ where r is the 
compression factor, which is a function of the Mach number ( r = 4M^/ (M^ — 1). 
Thus for a strong shock, (M >> 1), r = 4 and f(p) is proportional to , which 
corresponds to an energy spectrum proportional to E^^ at high energy. The 
observed proton spectrum {E~'^-'^) is slightly steeper due to an easier escape 
of high energy nuclei (see figure 5). Once launched in space, these nuclei are 
deviated by the magnetic field irregularities and they lose memory of their birth 
place. The propagation is of the diffusive type. Each time they reach the 
border of the Galaxy (badly defined) they can escape or return to the galactic 
disk where they suffer ionization energy losses and nuclear interactions. The 
Galaxy, in this context is like a leaky box. The influencial parameter is the 
mean escape length, A (or the related confinement time). At every moment on 
their way, the fast nuclei lose energy (by ionization and heating of the ISM) and 
can be destroyed in flight by a collision with other nuclei sitting in the ISM. 

The propagation equation reads 

dN(E)/dt = Q(E) - N(E)/t. d/dE [(b(E).N(E)] = 

(assuming a steady state). Q(E) is the source (injection) spectrum, b(E) is 
the energy loss rate by collisions with ambient electrons. The loss time r is the 
harmonic mean of the confinement time and the nuclear destruction time. The 
energy distribution is modified in the course of the propagation since energy 
losses are energy dependent. Dividing by the mean density of the ISM one get 

d<^(E)/dx = q(E) - 0(E)/A + d/dE[uj{E).cj){E)] = 

with uj{E) = dE/dx{MeVl{g.cm~'^)). 

The equilibrium spectrum (or interstellar spectrum) is solution of this equa- 
tion. It is over this spectrum and not the source spectrum that the cross sections 
have to be averaged to get the production rates and ratios. 

4.2 Low energy cosmic rays and the carrot 

At low energy (less than 100 MeV/n) the particles are thcrmalized before leaking 
from the Galaxy and being destroyed by nuclear collisions. Thus the propagation 
equation simplifies considerably since the Galaxy is now like a closed box. In 
this case, the production efficiency of LiBeB is maximal because fragments of 
nuclear collisions are quickly thcrmalized and integrated in the interstellar gas; 
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Figure 5: Spectrum of protons calculated through a diffusion-convection model 
[ |66| as a function of the kinetic energy. Solid line: model with no reacceleration 
and an injection spectrum index 2. Dashed line : interstellar spectrum of jTOf . 
Dash-dotted line: medium spectrum |Q , kindly supplied by Andrew Strong. 

on the contrary in the high energy case most of the LiBeB produced in flight 
escape from the Galaxy (only about 20 percent is remaining, ]6l|). 

The low energy part of the problem is often overlooked, whereas the key of 
the question could lie below 100 MeV/n since the peak of the cross sections are 
about at 20 - 30 MeV/n and even lower for alpha induced reactions. This is 
due to the large uncertainties affecting the interstellar flux of fast nuclei in this 
energy range. 

The propagation equation is now: 

q(E) + d/dE [w(E).^(E)] = 

Furthermore, as we said, the ^^B/^'^B ratio produced by GCR spallation, 
amounting to 2.5 (MAR, |6^) is somewhat different from the metcoritic value 
(4). MAR assumed the existence of a considerable excess of GCR at low energy, 
(5 - 40 MeV/n). This additional flux, called carrot, if existing, would enhance 
the production of ^^B via the reaction p + ^''N — > ^^B . 

It is worth noting that the composition of this carrot was taken as the GCR 
one and thus has to be considered as a secondary source and cannot be invoked 
to explain the linear evolution of B, indicated by the data. 

The difficulty is that this component cannot be observed in the earth vicinity. 
Indeed, below about IGeV/n, the energy distribution of GCR in interstellar 
space can only be determined indirectly since most particles are excluded from 
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the solar cavity by the solar wind and the theory of solar modulation is too 
uncertain to demodulate the observed flux, especially at the lowest energies. 

Limits on the low energy GCR flux (including the carrot) has been set 
by Lemoine et al. (1998) [ p9| running a galactic evolutionary model, in or- 
der to avoid overproduction of Be. A carrot of the kind assumed by MAR and 
Meneguzzi and Reeves (1975) is excluded since it would lead to strong over- 
production of LiBeB at solar birth. Only the lower spectrum proposed in the 
literature was consistent with the Solar System beryllium abundance. Apart 
this global (integrated) argument there are now rather secure means to get in- 
formations on the low energy spectrum. Strong and Moskalenko |66|, p3| 
have worked out a three dimensional model of GCR propagation and confine- 
ment taking into account all observables (cosmic-ray composition, gamma rays, 
electrons, positrons, synchrotron radiation and antiprotons,). They conclude 
that the low energy spectrum is depressed (figure 5). This conclusion has an 
impact on the production rate of the light elements since i) it excludes the carrot 
and ii) the decrease of the flux at low energy deduced confirms the restriction 
set by on the lower end of the spectrum and it leads to a good fit of the 
local LiBeB abundances. 

5 Are Galactic Cosmic Rays a truly secondary 
component? 

As shown above, the CR origin impacts on the galactic evolution of LiBcB 
elements. In the light of the linear relationship between BeB and Fe, certain 
authors have asked themselves if GCRs could be primary sources, at the expense 
of modifications in the physical picture of the origin of the accelerated nuclei. 

Indeed, since according to the reservoir of the accelerated nuclei (interstellar 
medium or supernova remnants loaded with fresh products of nucleosynthesis) , 
the production process of LiBeB is secondary or primary (i.e. dependent or 
independent of the interstellar metallicity Z). A pure ISM origin would lead to 
a slope 2 in the correlation between LibeB and O (or Fe) whereas the SN origin 
would lead naturally to slope 1 (see appendix, section 10). The debate on the 
origin and evolution of the light elements has shifted i) to the origin and very 
nature of the GCRs and ii) on the possible existence of a distinct low energy 
component (section 6). 

Concerning the status of GCRs, the question has recently been revived in 
the LiBeB Conference in Paris, December 1998 ]3^ . 

Maurice Shapiro representing the traditional trend, assumed that the 
cosmic rays are preaccelerated by coronal mass ejection (CME) driven shocks 
on low mass cool stars and accelerated further on by passing SN shock waves. 
Casse and Goret (1978) Q first pointed out that the difference between the 
compositions of the GCR source and the solar photosphere might be due to 
selection according to the first ionization potential (FIP) of the elements. The 
FIP effect is now well established as the cause for the differences between the 
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photospheric abundances and those of both the corona and the solar energetic 
particle in gradual events . The fact that these energetic particles are accel- 
erated by CME driven shocks in the corona is what drives the argument for the 
stellar origin of the GCR injection. However, arguments against this FIP based 
stellar origin were recently obtained from observations of ultraheavy nuclei in 
the cosmic rays p5| . 

According to Lingenfelter et al. (1998) |6^, the current epoch GCRs could 
come from each supernovae accelerating their own freshly produced refractory 
material. 

According to Meyer and collaborators, the current epoch GCRs originate 
from an average ISM of solar composition as in a classical proposal since the 
composition of the external layers of stars is identical to that of that of the ISM 
where they are born. Interstellar dust plays an important role in determining the 
abundances joj], |Q. In any case, an additional contribution from Wolf-Rayet 
circumstellar material enriched in freshly synthesized C and above all ^^Ne is 



needed to explain the Ne observed in GCRs |37 , [p3l. This WR component 



counts as a primary source, at least at present, but it is quite difhcult to assert 
its importance in the early Galaxy. 

The difficulties involved in the acceleration process of the SN model have 
been emphasized by Ellison and Meyer (1999) [Q. While the most significant 
shortcoming of this individual supernova model is the conflict with the delayed 
acceleration (see just below). Lingenfelter et al. (1998) show how the av- 
eraging of the supernova nucleosynthetic yields over the initial mass function 
and supernova types, as well as the inclusion of the effects of refractory dust 
grains formed in the ejecta, lead to a GCR source composition in good agree- 
ment with the observations. Moreover, Meyer and Ellison (1999) ||2^ stressed 
that the heavy s elements, and more specifically Sr, Zr, Ba and Ce, are not 
underabundant in GCR. They are thought to be produced in AGB stars and it 
is difficult to interpret this observation in terms of acceleration of SN material 
only. 

Observations of the electron capture radioisotope ^^Ni in GCR and its de- 
cay product ^^Co jl^, have been performed with an instrument on the ACE 
mission. ^^Ni decays by electron capture with a half life of 7.6 x 10^ years. 
However the decay is suppressed if the acceleration time scale is shorter than 
the lifetime because the atom is stripped as it is accelerated Q . The fact that 
much more ^^Co than ^^Ni is observed, suggests a delay (~10^ years) between 
explosive nucleosynthesis of iron-peak nuclei (under the form of their radioac- 
tive progenitors) and their acceleration to cosmic-ray energies. This makes it 
unlikely that individual supernovae accelerate their own ejecta. This conclusion 
is corroborated by energetic arguments. 

Obviously, the acceleration of average ISM matter is consistent with the 
implied delay between nucleosynthesis and acceleration. But this delay is also 
consistent with acceleration of ejecta matter in superbubbles where multiple su- 
pernova shocks accelerate accumulated supernova ejecta on average time scales 
at least as long as the implied delay ( pSj an d section 6) and are more energeti- 
cally favourable than single supernovae |47[ p8[. 
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6 Low energy nuclei from superbubbles: a pri- 
mary process 

As soon as the Hubble and Keck observations of BeB vs Fe were released, the 
necessity of a primary component has been felt |0 , , |^ . The existence of 
intense fluxes of fast C and O nuclei in the early Galaxy appeared unescapable. 
This LEG component is, in all likelihood, physically linked to the superbubble 
scenario. 

The physical conditions in superbubble cavities would lead to fast nuclei with 
hard energy spectra at low energies up to a cutoff at an energy which is still 
nonrelativistic due to a combination of weak reflected shocks and turbulence p2| , 
1^, H ||,|70[||. These are the low energy nuclei which have been postulated 



to produce the bulk of the Be and B at low metallicities On the other hand, 
as pointed out by Higdon et al. (1998) |5^, 0, since these giant superbubbles 
are thought to fill up a large fraction of the ISM, they are also the most likely site 
for the acceleration of the GGRs, which of course show no cutoff up to very high 
ultrarelativistic energies. Indeed the ingredients of a powerfuU ion accelerator 
are gathered in superbubbles due to copious injection of matter and energy 
by massive stars. WR and SN ejections in the common pot and supersonic 
material in the form of winds and supernova ejecta able to generate recurrent 
shock waves. Recent work by Parizot and Drury (1999) ||7^ shows that SBs 
produce Be and other light elements at the less energetic cost. They could be 
the key of the thorny energetic problem of the production of Be in the halo. 



Superbubbles |69 , | |55[ | thus appear as the most promising sites of acceleration 
of fresh products of nucleosynthesis, either as a separate low energy component 
or the GGRs themselves ||2^. Indeed, both the low energy nuclei and 

standard components could be produced by superbubbles at different stages of 
their evolution. 

This idea needs further observational substantation. Possible diagnostics 
of low energy particle interactions are i) non thermal X-rays in the 0.5 to 1 
keV range due to atomic deexcitation in fast O following electron capture and 
excitati on jsl, jsl [|| and ii) gamma- ray lines produced by nuclear excitation 
of G, O (72|, jeSl and '^Li -|- '^Be formed by the alpha -f alpha reaction ||^. For 
the time being, after the official withdrawal of the Orion gamma-ray results, |^ 
and the preliminary announcement of a detection of an excess between 3 and 7 
MeV in the direction of Vela |^ a hint to the existence of a large population 
of MeV/n nuclei comes from the X-ray emission of the galactic ridge ||8^. The 
observation or non observation of G, O lines at 4.4 and 6.1 MeV and of the Li-Be 



feature close to 500 keV by the INTEGRAL satelhte |125| wiU be the strongest 
test of the superbubble hypothesis |^ which are at the moment the best 
proposed sites of acceleration of low energy nuclei enriched in G and O. 
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7 Neutrino spallation : a primary process 



Neutrino spallation is a source of ^Li and ^^B via the interactions of neutrinos 
(predominantly j/^ and ly-r) on nuclei, specifically on ^He and ^^C [Q, Woosley 
and Weaver (1995) [0, hereafter WW. The lithium and boron yields are quite 
sensitive to the temperature of the fi and t neutrinos, in which there is a fair 
amount of uncertainty. As a result, the overall yields of ^Li and ^^B have 
considerable uncertainties, i^-process nucleosynthesis was incorporated into a 
model of galactic chemical evolution (Olive et al. (1994)^ and Vangioni- 
Flam et al. 1996 which have included the LEG component with the primary 
purpose of augmenting the low value for ^^B/^'^B produced by standard GCR 
nucleosynthesis. To correctly fit the observed ratio of 4, it was found that the 
yields of WW had to be tuned down by a factor of about 2 to 5 to avoid the 
overproduction of ^^B. Tuning down the z^-process yields ensured as well that 
the production of ''Li was insignificant, and did not affect the Spite plateau. For 
a review see p^ . 

Note that if taking the full yield, all galactic boron would be produced by 
p spallation. This could be a problem since ^Be is not coproduced and ''Li 
overproduced. Thus, this mechanism acts as a complement to nuclear spallative 
process at a level estimated to at most 20 percent concerning ^^B, if one wants 
to fulfill the observational constraints presented in section 2 pl[. 



8 Galactic evolution of light elements 

Analyzing all the physical parameters discussed above, two main LiBeB produc- 
ers emerge, the first one is the standard GCR in which fast p,a nuclei interact 
with CNO in the ISM. This process seems unable to produce sufficient amounts 
of LiBeB at the level observed in the halo stars. However, a recent study |pl[ 



[ p2| based on the O - Fe relation derived by | 111 ] and ^ at low metallicity (still 
controversial, see section 2) try to fit the observational constraints with a pure 
standard GCR (secondary production) component, but has difficulties with the 
B/Be ratio and possibly with energetic constraints at very low metallicity. 

The second one, LEG, invokes fragmentation of GO nuclei in flight by col- 
lision with H and He in the ISM. Massive stars are able in principle to furnish 
freshly synthesized G and O and accelerate them via the shock waves they induce 
in their surroundings. This mechanism is probably related to superbubbles. 

Finally, neutrino spallation is helpful to increase the ^^B/^'^B ratio up to the 
value observed in meteorites. It is also a primary process since it implies the 
break up of G within supernovae and not in the ISM. However it cannot be the 
unique mechanism to produce light elements since it does not make ^Be. 

These different mechanisms are included in a galactic evolutionary model to 
follow the whole evolution of each isotope. 

In brief, the characteristics of the standard galactic evolution model used 
[|[|0[ |l], |||l|, H are the following: 

No instant recycling approximation, in other terms account is taken of the 
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Figure 6: Beryllium and Boron evolution vs [Fe/H]. The halo evolution ([Fe/H] 
< -1, is dominated by the LEG component linked to massive stars. As far as B 
is concerned, there is room for a small contribution for v spallation. [M [E2| 



time delay of matter ejection by low mass stars. The cosmic ray flux and that 
of the LEG are taken proportional to the supernova rate (itself proportional to 
the star formation rate). The GGR production rate is time dependent through 
the growth of GNO in the ISM which is followed by the model (equiped with the 
relevant yields). The mass distribution of stars at birth, or initial mass function 
is proportional to M~^-^ from 0.4 to 100 Mo as usual. The star formation rate 
(SFR) is taken proportional to the gaz mass fraction all along the galactic life. 
The stellar lifetimes are taken from [116]. The composition of the ejecta is that 
calculated by WW. The model, of course, takes into account the destruction of 
LiBeB in stars. The SNI rate is taken constant, with the value observed today. 
This simple model is sufficient to map the evolution of each Light isotope in the 
abundance-abundance diagram (L-0 and L-Fe plots). 

The theoretical evolutionary curve of Be is normalized on the solar abun- 
dance of this element. 

Goncerning beryllium and boron, in this context, the main results are the 
following: the quasi-linearity (Be-B vs Fe) is easily reproduced thanks to the 
action of the LEG (fig 6, |93|). Standard GGRs contribute no more than about 
30 per cent to Solar System values. The B/Be ratio is in the range 10-30 as 
observed. The value 30 leaves enough room for neutrino spallation to reach 
iigyiiog _ ^ g^i- solar birth, but this contribution is marginal ]|9l|] . 

The ^Li/H ratio can also be explained in the framework of the same model 
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[ p2| , this without piercing the Spite plateau (figure 7). In this figure, showing the 
evolution of ^Li/H vs [Fe/H], it can be seen that GCR is overwhelmed by LEG. 
The decrease of the ^Li/^Be ratio could be explained in terms of the variation 
of the composition of superbubbles in the course of the galactic evolution, being 
O rich at start due to SNII and becoming more and more C rich due to the 
increasing contribution of mass loosing stars. 

Fields and Olive | ]5l| , |5^ using the new 0-Fe correlation and considering only 
the GCR component have also reproduced both solar and popll ^Li in a quite 
natural way. The proportions between the different production processes could 
be modified if the new 0-Fe relation is verified but even in this case, a primary 
component oparating in the very early galactic phase cannot be excluded. 




Figure 7: ^Li and ''Li evolution vs [Fe/H]. Horizontal full line: Li abundance, 
essentially ^Li. Full oblique line: LEG contribution to ^Li. Data points concern 
the two observed stars, (total Li and ^Li respectively). Dash-dotted line: GGR 
contribution to the ®Li abundance taking the classical flat 0/Fe relation. . 

In both scenarios, the evolutionary curve of ^Li crosses the halo observations 
meaning that ^Li is almost essentially intact in the envelope of stars in which it is 
measured. ''Li in turn, more tightly bound than ^Li, is even less destroyed, thus 
the mean value of the Spite plateau reflects nicely the Big Bang '^Li abundance. 
This reinforce the use of '^Li as a cosmological baryometer as shown also recently 



by the observational and theoretical analyses of |2S[] |100[ |104| 
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9 Conclusion and perspectives 



The confrontation of theory with observations, and the crossing of the bound- 
aries of experiment and theory, has permitted theorists, observers and experi- 
mentalists to clarify fundamental questions on the origin and evolution of light 
elements. Different evolutionary scenarios concerning LiBeB production by ac- 
celerated particles in the Galaxy, remain viable. 

1. The delay (^^10"'' years) between nucleosynthesis and acceleration, implied 
by the ^^Co and ^^Ni cosmic-ray data, is consistent with both acceleration of 
average ISM matter or ejecta matter in superbubbles. 

2. A low energy component distinct from standard GCRs appears necessary, 
specifically in the early Galaxy, but needs to be more observationally rooted. 

3. Future X-ray and gamma-ray line observations should help to clarify the 
very existence of the low-energy nuclei which are able to generate LiBeB in a 
primary mode, and if these low energy cosmic rays indeed exist, determine their 
energy spectrum. In spite of the problems encountered in the observation and 
analysis of the Orion region, the Vela region should be a prime target for the 
INTEGRAL mission. 

4. The new [0/Fe] vs. [Fe/H] observations may make it necessary to re- 
fine the model in which GCR nucleosynthesis, taken as a secondary process, 
is combined with nucleosynthesis by low energy nuclei originating from super- 
bubbles, which is a primary process. The contribution of each component will 
have to be adequately rescaled to account for the enhanced contribution of the 
secondary process; the contribution of neutrino spallation would also depend on 
this rescaling, but it will stay, in any case, marginal. 

5. New abundance measurements of ^Li close to [Fe/H] = —2.3, integrated in 
the framework of galactic evolutionary models, indicate that this nucleus is not 
significantly depleted in the observed stars and therefore Li (mostly ''Li which 
is less fragile than ^Li) is essentially intact in halo stars (Spite plateau). This 
makes Li an excellent indicator of the baryonic; density of the Universe. The ^Li 
Galactic evolution can be made consistent with the other light nuclei. 

Thus, the simplistic vision of the origin and evolution of LiBeB has devel- 
oped into a complex array of possibilities due to the wealth of observational 
discoveries. Nevertheless, the originally proposed fundamental mechanism re- 
mains preserved, namely the nuclear spallative process which is the only one to 
coproduce the light isotopes in the right proportion. 

The needs for the future are the following: 

On the theoretical side it would be necessary: 

i) to check NLTE corrections on B abundances (and perhaps Fe) at very low 

Z. 

ii) to develop and refine SN II models, especially at very low Z and high 
mass. 

On the observational side, it would be desirable to get measurements of ^Li, 
^Li, ^Be, i°B, "B, O, Fe in the same halo stars and to get ^Li/^Li and ^^B/^°B 
ratios in stars of both populations (popi and pop II). 
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Finally, the observation of C,0 and Li-Be gamma-ray lines are important ob- 
jectives of the INTEGRAL sateUite following the road open by the COMPTON 
GRO satellite. 

10 Appendix 

A clear distinction should be made between a primary production mechanism 

and a secondary one, in the sense of the galactic evolution. Indeed, in the 
plausible hypothesis that the GCR flux (f) is proportional to the supernova rate 
dN(SN)/dt, the rate of increase of L/H, L standing for beryllium or boron is: 

d{L/H)/dt = z{t){a<j){t)) 

where z(t) is the evolving CNO fraction by number and {a.<p(t)) the energy 
averaged of the production cross section times the flux. Since z is proportional 
to N(SN), the integrated number of SN up to time t, and (j) is proportional to 
dN(SN)/dt, we have, assuming a constant spectral shape, 

d{L/H)/dtaN{SN)dN{SN)/dt or L/Haz'^ 

On the contrary, freshly synthesized C and O nuclei injected/accelerated in 
SN ejecta or superbubbles would lead to a primary production of LiBeB, inde- 
pendent of the interstellar medium metallicity, by fragmentation of these species 
on the surrounding H and He nuclei. So the production rate is independent of 
metallicity z, and the cumulated abundances L/H are proportional to z. 

For example, to illustrate our purpose, we calculate the production rate of 
Be by protons fragmenting interstellar CNO by the GCR process: 

dN{Be)/dt = N{CNO)<j<t) 

Dividing by N(H), the hydrogen number density, one gets the rate of increase 
of the Bc/H ratio as a function of z defined as N(CNO)/N(H): 
d/dt(Be/H) = z{t).a(p{t) 

The present rate of increase of Be/H, calculated from realistic numerical 

values, N{CNO)/N{H) = lQ-^,a = 5m,h,(j) = lOcmr'^, for E > 30 MeV is 
about 10~^*s~^, which integrated over the galactic hfetime (10^° yrs) leads to 
Be/H about 10~^^ in qualitative agreement with the Solar System value, which 
is encouraging. However, this rough calculation does not take into account the 
temporal variation of the parameters. 
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